Spontaneous spin-polarized current in a nonuniform Rashba interaction system 
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We investigate the electron transport through a two-dimensional semiconductor with a nonuniform 
Rashba spin-orbit interaction. Due to the combination of the coherence effect and the Rashba 
interaction, a spontaneous spin-polarized current emerges in the absence of any magnetic material 
and magnetic field. For a two-terminal device, only the local current contains polarization; however, 
with a four-terminal setup, a polarized total current is produced. This phenomenon may offer a 
novel way for generating a spin-polarized current, replacing the traditional spin-injection method. 
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How to generate the spin-polarized current in a semi- 
conductor (SC) has been one of the most significant and 
challenging issues in condensed matter physicsi^ Apart 
from the fundamental physics interest, it may also have 
direct commercial applications. Over the past several 
years, the issue has attracted great experimental and the- 
oretical efforts. Due to the fact that semiconductors are 
in general spin-unpolarized, the key for generating polar- 
ized current in previous works is through spin injection, 
namely, to produce spin-polarized electrons from a po- 
larized source [e.g. Ferromagnet (FM) or polarized pho- 
ton], then inject them into a SC. However, among the 
currently existing spin injection methods^*^ none is very 
satisfactory. For the spin injection from a FM to a SC, its 
spin-polarization efficiency is usually low with a typical 
polarization around 1%& For the polarized optical meth- 
ods of spin injection, it is difficult for the integration with 
electronic devices £ 

Very recently, based on the Rashba spin-orbit (SO) 
interaction, some theoretical works have proposed dif- 
ferent approaches to generate a spin-polarized current 
without FM materialsi&LS However their devices are 
usually complicated or an external magnetic field is re- 
quired. The Rashba SO interaction is an intrinsic interac- 
tion in a two-dimensional electron system (2DES) of SC 
heterostructuresAifi It originates from an asymmetrical 
interface electric field, i.e. the asymmetrical potential 
energy in the direction perpendicular to the interface. 
The strength of the Rashba interaction can be tuned and 
controlled by an external electric-field or gate voltageii 

In this paper, we predict that a spin-polarized cur- 
rent spontaneously emerges in the SC in the presence 
of a nonuniform Rashba SO interaction. In particu- 
lar, this spin-polarized current is an intrinsic property 
of the nonuniform Rashba's SC and it does not need any 
magnetic materials nor a magnetic field. While under a 
voltage bias, a local polarized current is produced every- 
where, but with zero total polarized current. However, if 
for an open multi-terminal setup, a total polarized cur- 
rent emerges. Thus, our proposal offers an efficient and 
simple method to generate the spin-polarized current. 



We first show the principle of generating a spin- 
polarized current. For simplicity, we assume two paths 
for an electron traveling from one terminal of the sample 
to the other (see Fig. la), and t\ and t 2 are their respec- 
tive transmission coefficients. Because the Rashba inter- 
action strength a is tunable in experiments^ we choose 
different a in the two paths. An extreme case is a = 
in one path, e.g., path-1, and a large a in the path-2. 
This particular choice is not essential, but it brings out 
the physics more clearly. Due to the Rashba interaction, 
an extra phase is generated when an electron passes the 
path-2 In particular, this phase is dependent on the 
spin of the incident electron. For a spin-up electron, the 
extra phase is cp = —kuL = —am*L/h 2 (where L is 
the length of the path-2 and m* is the electron effective 
mass), assuming the Rashba energy is weak compared 
with the kinetic energy. On the other hand, the phase 
is — ip = kjiL for a spin-down electron. If only to con- 
sider the first-order tunneling process, the total trans- 
mission probability for the spin-up incident electron is 
Tj = \ti + t 2 e llp \ 2 , which in general is different from that 
for the spin-down electron, Tj = \t\ + t 2 e~ lip \ 2 . There- 
fore, a spin-polarized current is spontaneously generated 
and its polarization p at zero temperature is: 



2 l^iia j sin 8 sin tp 



?T +Ti \h\ 2 + \t 2 \ 2 + 2|tit 2 |cos0cos<p' 



(1) 



where 9 is the phase difference between t\ and t%. 

Next, we consider a specific two-dimensional and two- 
terminal SC system, shown in Fig.lc. In this device, 
two wires, II and III, are in the center region. In order 
to show our results are general, we choose the system 
without the mirror symmetry. In this set-up, an incident 
electron from Terminal I traveling to Terminal IV has 
two paths, i.e. passing the region II or III. If the Rashba 
interaction a's are different in the regions II and III, the 
above-mentioned coherent effect will occur. Then a spin- 
polarized current should be generated, although there is 
no magnetic material nor a magnetic field. 

The Hamiltonian for the two-terminal system (Fig.lc) 
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where V(x, z) is the potential energy. Here we let 
V(x, z) — in the region I and IV; V(x, z) = V2 (or 
V3) in the region II (or III); and V(x, z) = 00 in other 
regions. The last term in Eq.[2] is the Rashba interac- 
tion and a{x, z) describes its strength. For simplicity, we 
assume that a = in the regions I and IV, and a = 0,2 
and as in the regions II and III. Boundary matching is 
employed to solve for the transmission coefficients 
Assuming that the incident electron is at the subband n 
with the spin index s and the energy E from Terminal 
I, and to neglect the mixing of the inter-subband in the 
regions II and IlljlStli the wave functions §(x, z) in the 
regions I to IV are written as follows^ 



<&(x, z) = 

e ikI r> x (p^(z)s + r„ ms e-' lk ^ x ip I m (z)s 

m 
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s =| I I is the spin index, and s also describes the corre- 
sponding spin states, in which s = (1,0) T for f and s = 
(0, 1) T for |. yP m (z) {(3 = I, II, III, and IV) is orthonormal 
transverse wave functions for the subband m in the region 
/?. kii IV and k^ (7 = II or III) are the corresponding x- 

i/iv _ 



direction wave vectors with k, 
and feX„ 



± 



Vy - E2 



(e-eU iv ) 

sk Rj , in 



which k Rl 



a-,m 



</tf and Ei = ^(#£) 2 . t mns 



and r mns are the transmission and reflection amplitudes; 
a mns an< i ^mns are constants to be determined by match- 
ing the boundary conditions. Here the boundary con- 
ditions are i 1 ^ 1 ? <&(x, z)\ x=0 -/ L - = $(x, z)\ x=0 +/ L + and 
v x $(x, z)\ x=0 -/ L - = v x $(x, z)\ x=0+/L + + ^-$(Q/i, z), 
where v x = (p x + a z hk R )/m* is the velocity opera- 
tor and Uq is the Schottky 5 barrier potentials at the 
interfaces. 1 ^* Using the above boundary conditions, t mns 
can be exactly obtained, including all orders of reflec- 
tion and tunnelling processes. After solving t mns , the 
transmission probability T s can be obtained through the 



relation T S (E) = E m ,„ 9[E - E^E - E%) 
Similarly, the current (or conductance) density at an ar- 
bitrary location (x, z) can also be obtained. For instance, 
the conductance density gx s (x, z) in the x-direction in 
the region IV is: 



g Xs (x,z) 



dj 



Xs 



dV 



= — dE 



- df(E) 
dE 



where f(E) = \/[exp {E ~ EF V kBT + 1] is the Fermi distri- 
bution function with Ep being the Fermi energy. 
(2) We numerically study the conductance density 

gxs( x i z ) an d the local spin polarization p(x, z) = [gxi — 
9xi\/[gx] + 9xi\- I n the numerical calculations, we 
choose the system sizes to be: Wl = W R = L = lOOnm, 
W s d = Wn = 30nm, = 0, a R — SOnm, and 
D = lOnm. We also set k R 3 = and k R 2 = 0.015/nm, 
with the corresponding 0:2 = h J^* 2 ~ 3 x 10 _11 eVm for 
rn* = 0.036m e . Fig. 2 shows p(x, z) in the region IV. 
Here p(x, z) is clearly non-zero, and it can be over 15% 
at some locations. This means that the coherent effect as 
shown in Fig. la indeed plays the role for a finite p(x, z). 
For a further verification, we also study the following two 
cases for which the coherent effect is expected to vanish: 
(i) closing one channel, e.g. to make the region III to 
be very narrow; (ii) to set the a to be equal in both re- 
gions II and III (i.e. to set q?2 = a 3 ). Indeed, we find 
p(x, z) — in both cases for any (x,z). 

Now we show the behavior of the local spin polariza- 
tion p(x, z) in detail by plotting p(x, z) (the red dot- 
ted curve in Fig. 3c) and the corresponding conductance 
density gx]/i ( see Fig. 3a) versus z at x = lOOnm, i.e. 
the dotted line position in Fig. 2. Fig. 3a exhibits that 
gxT an d 9xi have clear difference. In particularly, at 
the peak position of gx\ n this difference remains, and 
it even reaches the largest value. Moreover the total 
conductance G s = J dzgx s (x, z) is quite large, ( e.g. 
Gj = G| w 1.2e 2 /h for the parameters of Fig. 2 and 
Fig. 3a, b). This means that this system can generate a 
large current density with a large local spin polarization. 
More importantly, the above property always survives, so 
long as the system size is within the coherent length. For 
example, at x = lOOOnm, gx] and gxi still have large 
difference (see Fig. 3b) and p(x, z) can exceed ±10% in a 
wide range of z (see the red dotted curve in Fig. 3d). 

Next we investigate how the local spin polarization 
p(x, z) depends on sample parameters, (i) When the po- 
tential V3 varies slightly, p(x, z) changes substantially. It 
can vary from the largest positive value to the largest 
negative value and vice versa (see Fig. 3c). This charac- 
teristic is very useful. Because V3 can be controlled by a 
gate voltage, so p{x, z) can also be tuned and controlled 
in an experiment, (ii) If there exists an interface poten- 
tial Uq, p(x, z) is barely affected. It may still exceed 10% 
(see the blue dash-dotted curve in Fig. 3d). But the con- 
ductances G s and gx s are weakened by a large Uq. (iii) 
With an increased distance D between the two channels, 
the overlap of two outgoing waves from the two channels 
is smaller, so p(x, z) will reduce slightly (see the black 
solid curve in Fig. 3d). But \p(x, z)\ can still exceed 5% 
for D = 50nm. (iv) With a larger Fermi energy Ep, 
more subbands in the region I-IV are available that in- 
creases G s and gxs- Meanwhile the variation of p(x,z) 
versus z exhibits a stronger oscillation^ and its ampli- 
s^l^ti^^e¥^fe4 v ma^gS , 'n1^ lashed curve in 



Fig.3d). 

(3) We emphasize although the local spin polarization 
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p(x, z) is fairly large almost everywhere^ the total con- 
ductance G s is unpolarization (i.e. Gf = Gj) for any 
two-terminal devices, because the two-terminal AB setup 
has the phase-locking effecti^S We prove the above state- 
ment in detail below. Due to the current conservation 
and the time-reversal invariance, the transmission coeffi- 
cient for a two-terminal AB system without the spin de- 
grees of freedom has the property of T(E, 4>) = T(E, — 0), 
the so called " phase-locking effect" , where cj> is the mag- 
netic flux through the AB loopS In our system, since 
no spin-flip processji 2 ^ the spin-up and spin-down elec- 
trons can be treated as two independent sub-systems. In 
the spin-up system, when an electron passes the lower 
channel, an extra phase tp — —kpi%L is added because 
of the Rashba interaction^ This extra phase plays the 
same role as if an external magnetic flux thread the AB 
loop. Then we have T^{E) = T(E, tp). Similarly, for the 
spin-down system, a fictitious magnetic flux —ip appears, 
and T i = T(E,-<p). Therefore, T r (E) = T^E) and 
G s = ^- J dE ~g£^ T s (E) must be spin-unpolarization, 
i.e. Gf = Gj for any two-terminal devices. 

In order to obtain a polarized total conductance (or 
current) , we devote the rest of the paper to study four- 
terminal devices. Consider a specific four-terminal device 
as shown in Fig. Id, in which the right (outgoing) termi- 
nal (the original region IV) is splitted into three terminals 
at the position x = L + hi. Assuming an incident elec- 
tron from Terminal I, the wave function $(a;, z) in the 
region I- VII (see Fig. Id) can be written similarly as for 
the two-terminal case. By matching boundary conditions 
at x = 0, L, and L + L^, the transmission amplitudes 
4ms( E ) ifi = V, VI, and VII) from the nth subband of 
Terminal I to the mth subband of Terminal [3 can be ex- 
actly obtained, although the deductive process is more 
complicated here. Afterwards, the transmission prob- 

ability TP{E) = £ ro ,„0(£ - E^E - O^flCJ 2 

and the conductance Gf = ^ / dE ^ 9 ^ E) T^(E) can also 
be calculated. In the numerical calculations, we choose 
the device geometry as: The left side and the center 
region II and III are the same as for the two-terminal 
device, and the sizes on the right are = 50nm, 
W 5 = W 7 = 50nm, W 6 = 30nm, Wr = 200nm, and 
L2 = lOOnm (see Fig. Id). To simplify, we set the poten- 
tial energy V and the Rashba interaction a in the region 
I, IV, V, VI, and VII to be zero. In a multi-terminal 
device, the total conductance Gf and the total current 
are spin polarized, so we focus on Gf and its polariza- 
tion pP [pi 3 = (Gf -Gf )/(Gf + G\ )}, instead of the local 
conductance gx s (x,z) and the local polarization p(x, z) 
as in the two-terminal case. 

Fig. 4a and 4b show the conductance Gjf 1 and its po- 
larization p VI versus the potential V3. G^ 1 and G^ 1 
show a large difference. This difference can be more 
than 0.15e 2 //i and p VI can exceed ±10% in a wide range 
of V3. p VI versus V3 exhibits an oscillatory behavior. 
In particular, it can oscillate from a maximum positive 
(or negative) value to a maximum negative (or positive) 



value with changing V3 . This characteristic is very useful, 
meaning that the spin-polarized direction and strength 
can be conveniently controlled in an experiment by tun- 
ing the potential V3. For the other two terminals V and 
VII, gY^ VI1 and p v ^ vn have similar behaviors. Below 
we emphasize two points: (i) It is the total conductance 
(or current) that is polarized, not only the conductance 
density with local polarization. This polarization can 
survive within the spin coherent length instead of the 
electron coherent length as in the two-terminal case. Usu- 
ally, the former may be much longer than the latter (ii) 
In the present device, the spin-polarized current is gener- 
ated without any magnetic material or a magnetic field. 
In the zero bias case, anywhere inside the sample is non- 
magnetic. When a bias is added, a spin-polarized current 
spontaneously emerges due to the coherent effect and a 
nonuniform Rashba interaction. 

We now study how the polarization p VI depends on 
other parameters: (i) p VI versus the Fermi energy Ep 
exhibits disorder-like oscillating behavior, and the am- 
plitude slightly weakens at high Ep (see Fig. 4c). (ii) p VI 
versus kp^ (i-e., a.2) is sinusoid-like curve with the pe- 
riod ~ 2ir (Fig.4d). But it is not exact periodic func- 
tion because the Rashba interaction also gives rise to 
an energy term h 2 k 2 R /2m* except for the extra phase 
— crfcijL. (hi) Fig.4e shows p VI versus has (i.e., as). 
Clearly £^3 = is not essential for a non-zero p VI . As 
long as \kn3 — kj^ \ ^ 0, a spin-polarized current appears. 

Finally, let us discuss the realizability. To add a gate 
(the deep gray region in Fig. lb) in a SC 2DGS, one 
can make the Rashba interaction a in this region dif- 
ferent from the a's in other regionsi^i Then under a 
bias, a local spin-polarized current is automatically in- 
duced. If four extra split gates (the black one in Fig. lb) 
are added to form an open multi-terminal device, a to- 
tal spin-polarized current is generated from source to 
drain. Notice that the device of Fig. lb has been real- 
ized about 15 years agoi 2 ^ Moreover, this device is much 
more open than the above-mentioned four-terminal de- 
vice (Fig. Id). The phase- locking effect is more severely 
destroyed, hence, this kind of set-up will have a much 
larger p. In fact, if the system is sufficiently open, then 
only the first-order tunneling process exists due to the 
current bypass effect, the spin-polarization p can reach 
100% at = \t 2 \ and 6 = tp = tt/2 [see Eq.(l)]. 

In summary, we propose a new method for generat- 
ing the spin-polarized current, replacing the traditional 
spin-injection approach. Here the spin-polarized current 
is induced due to the combination of the quantum co- 
herent effect and the Rashba spin-orbit interaction. In 
the two-terminal device, a local spin-polarized current is 
produced. While in an open multi-terminal setup, a to- 
tal spin-polarized current emerges in the absence of any 
magnetic material or an external magnetic field. 
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FIG. 1: (Color online) (a) Schematic diagram for an electron 
transport through two paths in a two-terminal device, (b) 
Schematic diagram for an open multi-terminal device made 
of semiconductor 2DGS with four split gates (the black re- 
gion) . The Rashba interaction in the deep gray region differs 
from the rest of the system, (c) and (d) are the configura- 
tions for the specific two-terminal and four-terminal systems, 
respectively. 

FIG. 2: (Color online) The local spin polarization p(x, z) 
vs x, z in the region IV for the two-terminal device. The 
parameters are V 2 = 0, V 3 = -0.02eF, E F = 0.013eV, 
U = 222^ = o, and k B T = 0. 

FIG. 3: (Color online) (a) and (b), g X \ (solid) and gxi (dot- 
ted) vs. z for x — lOOnm (in a) and x — lOOOnm (in b). (c) 
p(x, z) vs z for V3 = —0.023, — 0.02(the red dotted curve), 
-0.026, -0.028, and -0.031eV along the arrow direction, (d) 
p(x, z) vs z for the cases of (i) x = lOOOnm (red dotted curve); 
(ii) U — 0.2/nm (blue dash-dotted curve); (iii) D — 20nm 
and an = 20nm (black solid curve); (iv) E F = 0.05eV (ma- 
genta dashed curve). The others no-mentioned parameters 
in (a), (b), (c), and (d) are the same as for Fig.2 and at 
x = lOOnm. 

FIG. 4: (Color online) (a) Gj 1 (solid) and G\' (dotted) vs V 3 . 

(b) p VI vs V z for E F = 0.013eF (solid) and 0.015eF (dotted). 

(c) p VI vs Ef- (d) p VI vs k R2 for E F = 0.013eF (solid) and 
0.015eF (dotted), (e) p VI vs k m for k m = 0.015/nm (solid) 
and kn2 — 0.03/nm (dotted). The other non-mentioned 
parameters in (a)-(e) are V 2 = 0, V3 = — 0.02eF, k R2 = 
0.015/nm, k m = 0, and E F = 0.013eV\ 
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